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^T) ■ Abstract 

<N . 

We investigate the decoherence and dephasing of two coupled Josephson qubits. 
With the interaction between the qubits being generated by current-current corre- 
lations, two different situations in which the qubits are coupled to the same bath, 
or to two independent baths, are considered. Upon focussing on dissipation being 
caused by the fluctuations of voltage sources, the relaxation and dephasing rates 
are explicitly evaluated. Analytical and numerical results for the coupled qubits 
dynamics are provided. 

o 
o 

> 

^ ■ 1 Introduction 

Nano-electronic devices have recently attracted much interest as possible can- 
didates for the implementation of quantum computers. In fact, in contrast to 
ion traps [1] and NMR systems [2] , nano-electronic devices are easily embedded 
in an electronic circuit, and can be scaled up to large number of quantum bits 
(qubits) and quantum gates. Ultra small quantum dots with discrete levels, 
and in particular spin degrees of freedom embedded in nano-scale structured 
materials have been proposed [3]. However, these systems are difficult to fab- 
ricate in a controlled way. More suitable candidates are Josephson contact 
systems, where the coherence of the superconducting state can be exploited, 
and which can be fabricated by well established litographic methods. Macro- 
scopic flux states in a SQUID have been proposed as basic quantum bits [4]. 
Spectroscopy experiments of the flux qubit spectrum have recently been per- 
formed [5]. 
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In this work we consider the alternative design which uses the charge in low- 
capacitance junctions as quantum degree of freedom [6,7]. The logical states 
of the qubit are then macroscopic states differing by one Cooper pair charge. 
In a recent experiment Nakamura et al. have observed time-resolved coherent 
oscillations in such a Josephson junction set-up [8]. Several possible realiza- 
tions have been by now discussed [7,9,10]. The set up we consider is discussed 
in Section 2, and it follows the design proposed in [10]. In particular, single bit 
and two-bit operations may be performed by applying a sequence of suitably 
tailored gate voltages. 

Quantum computation requires a phase coherent time evolution. Hence, it is 
of prominent importance to deal with systems possessing long intrinsic phase 
coherence times, and to minimize possible external sources for dephasing [11]. 
Moreover, crucial for quantum computation is the ability to perform calcula- 
tions with single as well as with coupled qubits. Previous works analyzed only 
dephasing times of single qubits [12], or gave an estimate for N uncorrelated 
qubits [6]. Only a few previous works investigated the dynamics of coupled 
spins embedded in a thermal environment [13]. In this work we analyze the 
influence of dissipation and of thermal quantum fluctuations on the coherent 
evolution of two coupled qubits. We distinguish between two configurations: 
In the first one the two qubits are dephased by the same thermal bath. In the 
second one each qubit is coupled to its own environment. The general charac- 
teristics will be the theme of Section 3, where the phase coherence time and 
the decoherence rate will be explicitly evaluated for the two different environ- 
mental configurations. In Section 4 we provide analytical and numerical results 
for the population dynamics, when dissipation is dominated by fluctuations 
of the voltage sources. Finally, in Section 5 we draw some conclusions. 



2 Josephson Junction Qubits 

2.1 The Josephson- junction qubit 

In the following we briefly discuss the design for tunable Josephson charge 
qubits following [7,10]. It allows to manipulate the qubits parameters inde- 
pendently. The system consists of a tunable Cooper pair box, where a small 
superconducting island is connected to two tunnel junctions embedded in the 
two arms of a SQUID loop threaded by a flux $x, cf. Fig. 1. The junctions 
are supposed to be identical and are characterized by a capacitance Cj and 
coupling energy Ej. In this situation, an ideal voltage source V g is connected 
to the system via a gate capacitor C g . If the self-inductance of the SQUID 
is low, the two junctions effectively behave as a single junction being charac- 
terized by the effective capacitance Cj = 2Cj, and by the tunable Josephson 
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energy 



Ej(<b x ) = 2£°cos(7T$ x /$ ) , 



(1) 



where $o — h/2e is the flux quantum. Then, the charging energy of the 
superconducting island is characterized by the scale E c = e 2 /(2C g + 2Cj). In 
the following we consider the situation where the superconducting energy gap 
is the largest energy in the problem. Then only Cooper pairs tunnel through 
the superconducting junctions. In the regime Ej <C Eq a convenient basis 
is formed by the charge states parametrizing the number n of Cooper pairs 
on the island. In this situation the dimensionless gate charge n g = C g V g / (2e) 
acts as a control field. For values of n g close to a half integer two adjacent 
Cooper pair states become close to each other. We concentrate on the voltage 
interval near a degeneracy point where two charge states, e.g. n = and 
n = 1 play a role, while higher charge states, having much higher energy can 
be ignored. Hence, in the absence of dissipation, the system effectively reduces 
to a two-state system being described by the Hamiltonian 

H qh = ~E z (n g )a z - ^Ej($ x )a x , (2) 

where the cr's are Pauli spin matrices, and E z = AE c (2n g — 1) is controlled 
by the gate voltage. In this representation the charge states n = and n — 1 
correspond to the vectors | j) and | f) respectively. They are eigenstates of 
o z with eigenvalues —1 and 1, respectively. Hence, for a qubit being prepared 
in one of the eigenstates of a z , quantum coherent oscillations with frequency 
E q h/h, with 

£ qb := + Ej , (3) 

occur. 

2.2 The electromagnetic environment 

An ideal quantum system preserves quantum coherence, i.e., its time evolution 
is determined by deterministic reversible unitary transformations. In reality, 
any physical quantum system is subject to various disturbing factors which 
act in destroying the phase coherence. The concept of quantum computation 
heavily relies on the possibility of realizing quasi-ideal systems. Hence, it is 
crucial to investigate the effects of the environment on the qubits, and to 
understand how the environment induced dephasing can be minimized. It 
should also be realized that coherent quantum manipulations of the qubits 
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are still possible if the dephasing time is finite but not too short. In fact, the 
quantum error-correction techniques [15] allow to correct errors if they do not 
occur too often. 

In a charge Josephson qubit, the system is sensitive to the electromagnetic 
fluctuations in the external circuit and in the substrate, and to background 
charge fluctuations. We start by considering here the dissipative effects which 
arise from the fluctuations of the voltage sources. The equivalent circuit of 
a qubit coupled to an impedence Z{uj) is shown in Fig. 1. When Z(u>) is 
embedded into the circuit (with Ej = 0), the voltage fluctuations between the 
terminals of Z{lo) are characterized by the equilibrium correlation function 
[6,12] 



(5V(t)5V(0)) 

oo , 

= — J dujRe{Z t (uj)}JTW ( coth 



2k K T 



cos(c<jt) — % sin(cjt) 



(4) 



-l 



cj 1 



+ C n 1 , is the total im- 



Here Z t (u) := [iuC t + Z-\uj)\-\ with C t 
pedence between the terminals. Following [17] we model the dissipative influ- 
ence of Z t {uj) as resulting from a bath of harmonic oscillators described by the 
Hamiltonian 



#B = E 



pj 

2m, 



+ 



miU?x? 



(5) 



It is assumed that the voltage between the terminals of Z t {uj) is given by 

A 2 

5V = Y,i \ x i, an d the spectral density G{u) := | Y,i ~^ j zr&( u ~ u i) * s chosen 
to reproduce the fluctuations spectrum (4), i.e., G(uo) = ujReZ t (u). Thus, 
embedding the element Z(u) in the qubit circuit one arrives, in the two-state 
approximation, to the spin-boson Hamiltonian [6] 

#sb = # qb + H B + X -o z X , (6) 



where X := 2e(C t /Cj)5V describes the coupling to the bath. In the following 
we concentrate on the fluctuations due to an Ohmic resistor Z[uS) = R. It 
is convenient to introduce the spectrum Jy(c<;) := 4e 2 (Ct/Cj) 2 G(u)). For an 
Ohmic resistor it assumes the form Jv(oj) = 2nhauj which is linear at low 
frequencies up to some cut-off u c = (RCt)' 1 . The dimensionless parameter a 
which characterizes the strength of the dissipative effects reads 
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where Rk = h/e 2 pa 25.8kf2 is the quantum resistance. Hence, in order to make 
dissipation small one has to use a voltage source with very low resistance, and 
choose the gate capacitance C g ~ C t <C Cj as small as possible. For a typical 
resistance R pa 50f2 one has a pa 10~ 2 (C g /Cj) 2 . Upon choosing Cj pa IOOglF 
and C g pa laF one can reach coupling parameters as small as a pa 1CT 6 . 

A similar line of reasoning can be performed to include also the effect of 
1/ / noise due to background charges [16]. It is found that the effect of the 
background charges can be mapped into that of a harmonic oscillator bath 
with a spectrum that may have a simple behavior JyfiuS) oc K/uj. Usually, 
such fluctuations occur on a longer time scale than the voltage fluctuations. 
They also limit the coherence of the system. However in our charge qubit 
set-up dissipation is dominated by voltage flucutuations. In the following we 
shall mainly consider the effects of a harmonic bath being characterized by 
the Ohmic spectral density Jy(ou). However, how to include 1/f noise in our 
formalism will be discussed in the next Section . 

Finally, one can consider the effects of the fluctuations of the externally sup- 
plied flux through the SQUID loop of the charge qubit. These fluctuations 
couple to the qubit variable o~ x . We relate them to an effective impedence 
Ri, assumed to be real, of the current circuit. At typical high frequencies of 
the qubit's operation this resistance is of the order of the vacuum impedence 
Rj pa lOOfi. This yields for the dimensionless coupling constant [6] 



with M being the mutual inductance. For M ~ O.OlnH and Ej pa O.IK one 
obtains a x pa 10~ 8 . Thus, these latter fluctuations seem not to be dominant, 
and will be neglected in the following treatment. They can however be easily 
included in our formalism, as we shall discuss in Section 3. 

2.3 Coupled qubits 

For quantum computation pairs of qubits have to be coupled in a controlled 
way. In the scheme proposed in Fig. 2 all qubits are coupled by one mutual 
inductor L. For L = the system reduces to a series of uncoupled qubits, 
while they are coupled strongly when L — > oo. Hence, for finite L, the total 
Hamiltonian of this system consists of N qubits contributions as in (2), and 
of an oscillator resulting from the inductance and the total capacitance NCt, 
with C t _1 = Cj 1 + Cg 1 , of all qubits. The voltage oscillations in the LC t 
circuit affect all the qubits equally. They induce a phase shift 2-KLp = 2%^-^ 
of the phase variables 0j conjugated to the Cooper pair number rii of the 
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i-th qubit. Here $ denotes the flux in the mutual inductor, and $o = h/2e is 
the flux quantum. In the following the parameters are chosen such that the 
oscillator remains in the ground state for all relevant operation frequencies, 
i.e., E c ,Ej <C huj L , where oo L = (NLCt)^ 1 ^ 2 is the characterisitic oscillator 
frequency. Moreover, we assume that fluctuations of the phase shift ip are 
small compared to unit, i.e., 

< $ o- Since ($ 2 )/L « hw L /2 this 
condition imposes only weak constraint on the parameters. Then, although 
the LC t oscillator remains in its ground state, it provides an effective coupling 
between the qubits. It reads [10] 

H^ = -Y, Ej{ * x f J ^ Xj) ^\ 0) 

i<j 



with the energy scale 




^=l7fl^. (10) 



while E) = Ej(<fr X i) are the effective Josephson energies of the qubits, con- 
trolled by the external flux. This coupling energy can be easily understood as 
the magnetic energy of the current in the inductor, where the current is the 
sum of the contributions from the qubits with non zero Josephson coupling, 
li oc Ej(j( l \ Noticeably the strength of the interaction does not depend on 
the number of qubits in the system. However, the frequency wl of the os- 
cillator scales with l/y/N. This limits the allowed number of qubits in the 
system, since this frequency should not drop below typical eigenfrequencies 
of the qubit. The interaction scale El involves the screening ratio C g /Cj. 
As discussed in the preceding subsection, this ratio should be taken as small 
as possible to minimize the decoherence effects of the dissipative electromag- 
netic environment. Consequently, to achieve a reasonable interaction strenght 
a large inductance is needed. For example, for the typical values Ej m lOOmK 
and Cg/Cj ~ 0.1 an inductance L m lfiB. is needed. A more sophisticated 
set-up which overcomes this problem enabling to use smaller inductances is 
discussed in [6]. 

From the above discussion it is clear that, with tunable Josephson couplings 
Ej(&x), two-qubits gate operations can be performed by setting to zero all the 
couplings except for two selected qubits, say 1 and 2. The resulting two-qubits 
Hamiltonian takes the form: 



Pl /?2 rpl rp2 rpl rn2 

TJ J J 2 ^(2) J J „(1) „(2) /-,-,>> 

# 2 qb = — 2~0* ~ ~Y°x ~ -y^ - -y^ E~ a v a v • ( U > 
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Finally, to include the environmental influence, we distinguish between the two 
different situations shown in Figs. 2a and 2b. In case I, shown in Fig. 2a, all 
the qubits are coupled to the same oscillator bath. The resulting Hamiltonian 
then reads: 

Hj = H 2qh + H B + + aW)X , (12) 

with X defined below (6). As we shall see in the next section, this case reveals 
interesting symmetries. In case II, shown in Fig. 2b, each qubit is coupled to 
its own oscillator bath. It is described by the Hamiltonian 

Hn = E #sb + hP*? , (13) 

i=l,2 Z 

with Hsb being the spin-boson Hamiltonian introduced in (6), and where the 
coupling parameter 

7 = -2^3 (14) 

was introduced. This case of uncorrelated baths is thought to better describe 
the situation in charge Josephson qubits. 



3 Dissipation and dephasing of coupled qubits 



In this section we explicitly solve the dynamical problem posed by the dissi- 
pative two qubit Hamiltonians (12) and (13). To better outline the symme- 
tries and possible constants of motion of the problem, we find it convenient 
to discuss the coupled qubit evolution in the Hilbert space spanned by the 
triplet states | TT) := (1,0,0,0) T , (| Tj) + | |T))A/2 := (0,1,0,0) T , together 
with | ||) := (0,0, 1,0) T , as well as by the singlet state (| ||) - | !T))A/2 : = 
(0, 0, 0, 1) T . In this representation the dissipative Hamiltonian for uncorrelated 
baths assumes the appealing form 



#TT = 



e — s 


V 


7 


-At] 


V 


-7 


V 


-Ae + As 


7 


V 


— e + s 


Ar] 


-Arj 


-Ae + As 


Arj 


7 



(15) 
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Here, for simplicity of notation, we introduced the quantities 



As := El - E 2 Z , Ar] := (E) - Ej)/ \pl , 
e:=El + El, r] := (Ej + Ej)/ \pl , (16) 

characterizing the isolated coupled qubits. The sum and difference 

s:=X 1 +X 2 , As = X 1 -X 2 , (17) 



characterize the baths coupling, cf. below (6) for the definition Xi for the 
single qubit %. 

For the Hamiltonian Hj of qubits coupled to a single bath the same form as 
for Hu applies upon substitution of 

s -> 2X , As = . (18) 



3.1 The degenerate case 



In the following we focus on the the degenerate case Arj = Ae = 0, corre- 
sponding to equal Josephson energies and asymmetries of the two qubits. It 
is in this situation that quantum effects arising from the coupling between 
the qubits are expected to play a major role. Then the nondissipative Hamil- 
tonian H 2q b becomes separated into two blocks, corresponding to the triplet 
and singlet states, respectively. To be definite, the unperturbed Hamiltonian 
satisfies the Schrodinger equation H 2q b\n) = E n \n), with eigenvalues 

£ M = ^^ 2 + f + 2f :=±|, £2,3 = ^, (19) 
and eigenfunctions, expressed in the singlet /triplet basis, 



|1) = ^ , 1 -^ (E + 1 + e,2r ] ,E-e + 1 ,0) T , 

2y/E(E + 7) 

|2> = (0, 0, 0, 1) T , |3) = ^=^ (V, -e, -V, 0) T , 
2 JE(E - 7) 



Hence, if the system was prepared in one of the triplet states, it will undergo 
quantum coherent oscillations among the three states with oscillation frequen- 
cies uj mn := \E m — E n \/h, with m, n — 1, 3, 4. However no motion is associated 
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with the singlet state, being the eigenstate |2) of the unperturbed Hamiltonian. 
We observe that the maximal oscillation frequency is U41 = E/h> V2E qh /h, 
where E^ is the single qubit tunneling splitting (3). In the case of distinct 
baths, dissipation mixes the triplet and singlet states. To be definite, the total 
Hamiltonian for nncorrelated baths reads H u = H u + H B with 



S 7] 



7 







i] —7 i] As 
7 i] —e + s 



V 



As 



7 ) 



(21) 



Because the problem of coupled qubits is isomorph to that of two interacting 
spins in the presence of dissipation, it is suggestive to associate to the triplet 
and singlet states quantum numbers (S, S z ), being eigenvalues of the squared 
total spin operator S 2 , and of the projection S z of the total spin on the z-axis. 
Then, the triplet states are characterized by \S — 1, S z — — 1, 0, 1), while the 
singlet state is given by \S — 0, S z — 0). Hence, once the system is prepared, 
e.g., in the entangled state with quantum numbers \S — 0, S z — 0), it will stay 
in that state for ever if no dissipation is present. In fact, as appearent from 
inspection of H^, the effect of dissipation is to induce transitions between the 
states \S — 0, S z — 0) and \S = 1,S Z = 0), yielding thermalization also for 
the singlet state. 

Thermalization of the singlet state however is impeded in the case of dissipa- 
tive qubits being characterized by the Hamiltonian Hi = + Hb- In fact, 
in the degenerate case is obtained from upon performing the substitu- 
tion (18). When dissipation is added, the oscillation frequencies u nm , m 7^ n, 
acquire a finite dephasing rate T mn . Moreover, the system relaxes to thermal 
equilibrium with a decay rate T being associated with incoherent tunneling 
processes. These rates are explicitly evaluated in the next section. 



3.2 Redfield equations 



To perform quantitative calculations we apply the well established Bloch- 
Redfield formalism [18,19] to the Hamiltonians H\ and Hu. The Redfield ap- 
proach provides a set of coupled master-equations for the matrix elements of 
the reduced density matrix. Upon performing a Markov approximation, and 
in the basis of the eigenstates of the unperturbed Hamiltonian, they read 

Prim i-^nmPnm ^ ' R-ra,m,n/ ",mi Pnlml j (^2) 
nl,ml 
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where R is the Redifield tensor whose elements are given by 



R-n,m,n/,ra/ ^-ml,m,n,nl + ^-nl,n,m,ml \^-n,k,k,nl^m,ml + ^- m ,k,k,mt^n,n/ ■ (23) 

k 

For the Hamiltonian if n the tensor A assumes the form 



v n,m,nl,ml ~ fc^> i S nm s n/ml + r nrn r ntm/) M\i(u) n i m i) 



~ ^{ s nm r nlml + r 'nm s nlm/) AM (uJ n , m/ ) , (24) 

where, M n (fi) := Mi (ft) + M 2 (ft), AM(fi) := Mi (ft) - M 2 (ft), with M;, 
(i = 1,2), reading 

oo 

Mi = -J dte- iQt (X t (t)X t (0)) . (25) 
o 



Upon introducing the spectral densities Juioj) := [Jvi{oj) + Jy 2 (u;)]/4 and 
AJ(u) := [J yi (cj) — Jy 2 (u;)]/4, cf. below (6) for the definition of Jy for the 
single qubit, the one-sided Fourier transforms can be explicitly evaluated from 
(4). One finds 



lmMn(n) = -7du,- Muj) 



7T 



h 



ReM„(ft) = -J„(ft) 




(26) 



where u c is the characteristic cut-off of the bath spectral density J n , cf. below 
Eq. (6). Likeways AM(w) has the same form as (26) but with J\i(u) — > AJ(w). 
In the following, for simplicity, we focus on the case of equal baths, yielding 
AM = 0. Finally, s mn and r nm are the matrix elements of the operators 



s = 



a z + °i )/2 and f = (a)P — af>)/2. We note that s nm derives from 
transitions within the triplet space, while r nm describes transitions between 



the singlet and the triplet state with the same quantum number S z 
be definite these matrix elements read: 



0. To 



Sll = -S 4 4 = , S 2 2 = S33 = , S13 = . = 



S34 



yfE&f + e 2 ) ' 



S14 



^^(2^ + £ 2 ) ' 



(27) 
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Snm = s mn and S2j = Sj2 = for any j. Finally, the only non vanishing matrix 



elements r nm = r rnn are 



ri2 = * , r 2i = , , r 23 = , (28) 



<E{E + i) y E(E — 7) y/frf + e 

For the Hamiltonian Hj the tensor A reads 



^n,m,nl,ml ~ ^ s nmS n imlMi(uj n , m ,) (29) 



where, 



oo 

Mi 



(Q)= J dte- int (X{t)X{0)) . (30) 
o 

Hence, Mj(f2) has the same form as (26) with Ju(u) — > Ji(cj) = Jy(uo). Note 
that, for the case of equal baths here considered, 

Ji(u) = 2J n (w) = 27rnac; , (31) 



with a defined in (7). Thus, apparently, dissipative effects are stronger for the 
dissipative configuration I. Some comments are necessary: 

(i) These results are valid for a dissipative coupling to the qubit variable o~ z . 
A generalization which includes dissipation parallel to the other spin matrices 
crj, as for example the flucutuating contribution (8), is straightforward. In fact, 
such contributions simply modify additively the tensor A njmjn/jm/ . It acquires 
extra terms which are proportional to the matrix elements of the operators 
o~\ ^ ± af^ , and of the spectral density of the new dissipative contribution. 

(ii) As discussed in the previous section, together with voltage fluctuations 
possessing an Ohmic spectrum, the electromagnetic environment provides a 
1/ / contribution arising from the background charges. This noise also couples 
to the a z operator of the single qubit. Hence, to take into account this noise, 
as well as other noise sources coupling to cr z , the total spectral density for the 
single qubit capturing the effects of a z coupling is modified to 

J z (co) = J v (uj) + J 1/f (w) + ... (32) 



This introduces different time scales for the decay of noise-induced correla- 
tions. In particular, those associated to 1/f noise decay on a very long scale 
compared to the qubit dynamics. Hence, the Markov approximation invoked 
to obtain the set of equations (22) breaks down. The generalization of (22) to 
include non-Markovian effects should be considered. 
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3.3 Relaxation and dephasing 



Due to the environmental influence, the oscillatory motion of the undamped 
system gets damped. Moreover, a frequency shift of the bare oscillation fre- 
quencies u) nm occurs. We evaluate the dephasing rates T nm and the decoherence 
rate T to lowest order in the coupling to the bath from the Redfield equations 
(22). We find the expressions (i = I, II) 



4 

I 

n>m=l 



r« = -]TR„, n ,„, n (A*) = 2 £ (33) 



with Si(u> nm ) := [Ji{u nm )/h\ coth(hw nm /2k B T), and 7 * m := s 2 nm + r 2 nm , 7 ^ m : = 
s 2 nm . Likeways for nondegenerate levels (|u; mn | > \R n ,m,n,m\) the dephasing 
rates read 



Tf m = -ReR ri , m , n , m (A*) = ^(s nn - s mm ) 2 M l (0) + 7 * tn 5 , i (w nm ) 



(34) 



Thus, the real part of the Redfield tensor provides the relaxation and dephas- 
ing rates. The imaginary part is responsible for an environment induced shift 
of the oscillation frequency u mn . One finds u nrn -> u mn := u nm - ImRn^^, 
where 



- [iuR, ,„ , ,„ = {s 2 nn - s 2 mm )lmMi{0) + l l nm ^- [ du ^\ coth ' 



■n,m,n,m ' 



q mn 



2k B T 



We observe that at zero temperature the general relation holds 



(35) 



T« = ^—r E r nL T = , (36) 



where iV; ei , is the number of levels involved in the dynamics (iVj e „ = 4 for cou- 
pled qubits, and N iev = 2 for uncoupled qubits). It follows that the dephasing 
rates are always smaller than the relaxation rates T^ l \ It is now interesting to 
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compare (33) and (34) to the decay rates T q b and describing relaxation and 
dephasing, respectively, of the uncoupled single qubits. An analysis similar to 
the one here reported yields for the single qubit the results [12] 

r qb = ^rS(E qh ) , r, = ^ + ^-5(0) , (37) 

where S(lu) = [Jy(uj)/h] coth(ftw /2k B T). As shown by (33), incoherent pro- 
cesses contributing to increase with the number of levels involved. Thus, 
the coupled qubits relaxation rates are larger than the single qubit relax- 
ation rate T q b, cf. also next Section and Fig. 4. As far as regards the dephasing 
rates, only the last term in (34) is not present for a single qubit. This term, 
however, is expected to play a minor role. Hence, the dephasing rates T mn 
remain of the same order of the single qubit relaxation rate T^, cf. Fig. 4. 



4 Dynamics of dissipative coupled charge qubits 

By gathering all of the results of the previous Sections 2 and 3, we evaluate 
now the occupation probabilities P^{t) := p w (t), with the vector G {| ft 
)> I II)) (I Tl) =t I IT))/V2}- Given an initial reduced density matrix p(t ), we 
obtain with t = and to linear order in the coupling parameter a, the result 

P,(t) = E «Jn[Pnn(io) - Pnn]e- rt + P? 
n 

4 

(38) 

where a Mn := (fi\m) are the matrix elements for the transformation from the 
eiegenvectors basis to the singlet /triplet representation, cf. (20). Moreover, 
the stationary equilibrium value is 

P? = E O/mWC • (39) 

n,m 

Note that, due to dissipation, the offdiagonal elements p^ m := p mn (t = oo), 
n 7^ m, are different from zero. They are terms of order a and vanish with 
a -> 0. 

In the following we present results of a numerical integration of the Redfield 
equations (22). We consider the case of identical qubits coupled to the same 
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bath with spectral density Jviyj) = lnhau). For the discussion of the results 
we find it convenient to introduce the frequency uj := Ej/h, where Ej = 
E] = E]. We start by showing in Fig. 3 a comparison between the analytical 
prediction in Eq. (38) and results of the numerical evaluation of the Redfield 
equations. In particular, we plot in Fig. 3 the time evolution of the survival 
probability P^(t) of finding the system in the triplet state | f f ) in which it was 
prepared at the initial time to — 0. The different oscillatory bits corresponding 
to the frequencies uj nm are clearly seen. A good agreement with the analytical 
solution Eq. (38) is found. We observe that, since in (38) only linear terms 
in a are kept, deviations from the numerical solution increase with increasing 
of the coupling parameter a (not shown). We also note that, to check the 
validity of the Markovian approximation involved in the derivation of the 
Eqs. (22), a numerical integration which abandoned the Markov assumption 
was performed. In the chosen parameter regime a perfect agreement with the 
Markovian- Redfield result was found (not shown). Figs. 4 shows results for 
the survival probability P^(t), and for the occupation probability P e nt(0 °f 
the entangled state (| fl) + I !t))/\/2 for two different values of a. The 
initial condition p(t ) = | j"T) is assumed. The figure clearly shows a strong 
damping of the coherent motion already for the "small" damping constant 
value a = 0.01. However, when a = 0.001 coherence is preserved over many 
coherent oscillation periods. Note also the phase shift induced by the different 
choice of the damping parameter a. Fig. 5 shows the relaxation rate V (dashed 
line) and dephasing rate T 14 (solid line), versus temperature. For comparison 
also the (37) relative to a single qubit is shown. Clearly, the single qubit 
rates turns out to be smaller than its coupled qubit counterpart, but it is 
still of the same order of magnitude. When the qubit relaxation rate T q b is 
compared to the coupled qubits decay rate T, however, a larger increase of T 
in comparison to T qb is found, as expected from (33) (not shown). 

Finally, in Fig. 6 results for the relaxation and dephasing rates versus the cou- 
pling parameter —7 are shown. To better outline the differences with respect 
to the single qubit case, the rates have been normalized to the dephasing rate 
r^, of a single qubit. Note that in the absence of coupling is Ti 4 = Ti 3 = 2r^. 
Moreover, the relation (36) always holds. 



5 Conclusions and outlook 

In conclusion, we evaluated the decoherence and dephasing rates of coupled 
qubits for two possible dissipative configurations. The first case, denoted I, 
dealt with the case in which the coupled qubits experience the same dissi- 
pative forces, arising from a common electromagnetic environment. In the 
second case, denoted II, each qubit is coupled to its own oscillator bath. Due 
to enhanced correlations, dissipative effects are stronger in configuration I. 
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Interestingly enough, we found that in the degenerate case, in which the two 
qubits are characterized by the same Josephson coupling energy and the same 
asymmetry energy, the nondissipative Hamiltonian becomes separated into 
two blocks, corresponding to triplet and singlet states, respectively. Hence, a 
system prepared in one of the triplet states undergoes quantum coherent os- 
cillations among the three states, while no motion is associated to the singlet 
state. When the coupling to the correlated environment I is included, thermal- 
ization of the triplet states occurs, while the singlet state remains unaffected. 
In contrast, for the case of uncorrelated baths II, environmental induced tran- 
sitions between the singlet and triplet states occurs. Hence, also the singlet 
state thermalizes. 

When compared with the relaxation and dephasing rates of a single qubit, it 
turns out that the coupled qubits relaxation rate scales with the number of 
allowed interlevel transitions. This yields an enhancement of decoherence ef- 
fects as compared to the single qubit case. We performed explicit calculations 
which kept into account the effects of fluctuation of the voltage sources. In 
fact, this is thought to be the dominant decoherence mechanism in Josephson 
qubits. However, the effect of other noise sources can be esaily icluded in our 
formalism. 
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Fig. 1. Charge qubit with Josephson coupling controlled by an external field <&x 
threading the dc-SQUID. The dissipative influence originating by voltage fluctua- 
tions are captured by an impedence Z(u). 
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Fig. 2. A register of many charge qubits coupled to a common inductor L. In case a) 
the qubits are subject to a common dissipative environment, while in configuration 
b) each qubit is coupled to its own oscillator bath. 




Fig. 3. Time evolution of the survival probability P++ for coupled qubits being 
prepared in the state | jj). The solid line denote results of a numerical integration 
of the Redfield Eqs. (22), while the dotted line refer to the analytical solution (38). 
We choose Ohmic constant a = 0.0001 and zero temperature. Here and in the next 
figures we introduced ujj = Ej/h and we set E\ = E\ = 0, 7 = — 0.1-Ej, u> c = 50wj. 
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Fig. 4. a) Time evolution of the survival probability P^, and b) of the probability 
of the triplet entagled state P en t(t) for coupled qubits being prepared in the state 
ft). The curves, obtained from a numerical integration of the Redfield Eqs. (22), 
show results at T = for coupling parameter a = 0.001 (dashed line), and a = 0.01 
(solid line). 
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Fig. 5. Relaxation rate T and dephasing rate Ti4 of coupled qubits versus tempera- 
ture. For comparison also the single qubit dephasing rate is shown (dotted line). 
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Fig. 6. Relaxation rate T and dephasing rates r34,Ti3, and Ti4 of coupled qubits 
versus the qubits coupling parameter —7. The rates are normalized to the single 
qubit dephasing rate IV The temperature is set to zero. 
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